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We investigate the time evolution of the temperature and entropy of a gravitationally collapsing de 
Sitter Schwarzschild domain wall as seen by an asymptotic observer. Recent work has completed this 
analysis for a Schwarzschild domain wall and for a 3+1 BTZ domain wall. There were some striking 
qualitative differences between the two. Specifically, the BTZ domain wall exhibited a decrease in 
entropy over time, implying that it will not collapse spontaneously. However, the BTZ domain wall 
contained both a cosmological constant and a different topology from the Schwarzschild domain 
wall, and we wish to isolate which of these is responsible for the qualitative differences. Hence, we 
will study the de Sitter Schwarzschild domain wall, as it has identical topology to the Schwarzschild 
domain wall yet also contains a cosmological constant. We utilize a wavefunctional approach where 
we couple a scalar field to the background of the collapsing domain wall and determine the spectrum 
of the radiation as a function of time. The spectrum is quasi-thermal, and the degree of thermality 
increases as the domain wall approaches the horizon. The fact that the distribution is thermal allows 
for the determination of the temperature as a function of time. Since the collapsing domain wall 
forms a black hole, we compare the late time temperature to the Hawking temperature and find 
very good agreement. The temperature also exhibits the proper scaling with the mass of the domain 
wall. Using the thermodynamic definition of entropy and the temperature-mass relation, we find the 
entropy as a function of time. We find that the temperature profile exhibits an interesting feature 
in the beginning of its collapse where the temperature increases for a short time before decreasing 
and approaching a constant. The entropy exhibits a similar feature in that there is a very small 
decrease in entropy before it starts to increase and eventually approach a constant. These features 
go away for smaller values of the cosmological constant. This suggests that 1) gravitational collapse 
is prevented in a certain parameter space, and 2) the probable cause of the qualitative differences in 
the domain wall collapse of the Schwarzschild versus BTZ domain walls is the cosmological constant 
and not the topology. 



I. INTRODUCTION 



The properties of static black holes are currently very 
well known. This work was pioneered by Bekenstein, 
Gibbons, and Hawking. They showed that quantum pro- 
duction of particle-antiparticle pairs near the horizon of 
a black hole could make it appear as though the black 
hole were radiating. It was also determined that the 
spectrum of this radiation was thermal and that the en- 
tropy of the black hole was proportional to its area [T]- 
|3]. The conventional process to determine the entropy 
uses the Bogolyubov method, which works by determin- 
ing Bogolyubov coefficients that connect the initial vac- 
uum state far outside the matter to the final vacuum state 
after the black hole is formed. As is evident, this method 
only utilizes the initial and final states of the system, and 
therefore there is no knowledge of the time dependence. 
Recently, Ref.[lH3 developed a method for determining 
the spectrum of radiation given off by a collapsing do- 
main wall in a time dependent manner. These papers 
used a spherically symmetric collapsing domain wall for 
their analysis, and Ref. TIT applied it to a 3-1-1 BTZ do- 
main wall. The most striking qualitative difference be- 
tween the behavior of the two types of domain walls was 
that the entropy of the BTZ domain wall decreased over 
time, suggesting that the gravitational collapse would 
not occur spontaneously. The 3-1-1 BTZ metric, however, 
had two key differences from the Schwarzschild domain 



wall. The first and most obvious difference was a different 
topology. The second difference was that the BTZ do- 
main wall metric contained a cosmological constant. For 
this reason, we decided to determine the time evolution 
of the entropy of a de Sitter Schwarzschild domain wall, 
as it has identical topology to the Schwarzschild domain 
wall but contains a cosmological constant. To complete 
the analysis of the time evolution of the temperature and 
entropy, we will first determine the equation of motion 
for the collapsing domain wall from an asymptotic ob- 
server's point of view. This equation of motion will be 
used in the quantization of a wavefunctional that contains 
a scalar field coupled to the background of the domain 
wall. Then, using the t=0 wavefunctional as a basis, we 
will determine the occupation number as a function of 
frequency. This in turn will allow for the determination 
of the thermodynamic quantity /3 and therefore the tem- 
perature. We then determine the late time temperature 
as a function of horizon radius. Using the thermody- 
namic definition of entropy, dS = dQ/T, we integrate to 
find entropy as a function of temperature, and since we 
know the temperature as a function of time we can find 
the entropy as a function of time. We compare our results 
against the standard Hawking results, and we comment 
on the results. 
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II. EQUATION OF MOTION OF THE DOMAIN 
WALL 

It was shown in Ref.0] that for a domain wall with a 
general outside metric 

ds'' = -fir)dt^ + -^dr^ + r''dn'', (1) 

the classical equation of motion near the horizon from an 
asymptotic observer's point of view is 



i?«±/ 1 



2 



(2) 



The dynamics for R ^ Rh can be obtained by solving the 
equation R = ±f. For a de Sitter Schwarzschild metric, 



/=1- 



R. 



Ai?2 



(3) 



R 3 

and the equation of motion to first order in Rq — Rh is 

Rit) ^Rh + (i?o - i?H)e-(i/^«-^^-)*, (4) 

where Rh represents the black hole horizon radius. It 
should be noted that for a de Sitter Schwarzschild black 
hole, the horizon radius Rh is not Rg — 2GM. In fact, a 
de Sitter Schwarzschild black hole has two horizons. The 
outer horizon is called the cosmological horizon, and in 
order to produce the proper dynamics we must ensure 
that the collapsing domain wall starts within the cosmo- 
logical horizon. 

III. RADIATION WAVEFUNCTIONAL 

First, we will consider the radiation given off by the 
domain wall during gravitational collapse. This will be 
done by coupling a scalar field to the background of the 
collapsing shell. We will decompose the scalar field as 



$ = ^ak{t)ukir). 



(5) 



The exact form of Uk{r) will not be important to us. To 
find the modes ak{t), we will insert the metrics given by 
Eqs.(??) and (??) into the action 



(6) 



where we are using r — R{t) as the domain wall's posi- 
tion. The total action will be written as the sum 



S = S, 



in ~t~ Sout 



(7) 



where 



S,.n = 2tt I dtj^^^'^ drr^ (^-^^^ + T{dr<S>)^ ] (8) 



is the metric inside the shell, 
Sout = 2tt dt drr^ 

J JR{t) 



f 



+ fidr'^r (9) 



is the metric outside the shell, and 

dT _ I I [A-nR-^ ' 

According to Ref . [T^ . 



fR^ 



where 



h = 



so dT/dt can be rewritten as 
dT f 



(10) 



(11) 



(12) 



i?4 



dt AV^ + [^--^^/. 



(13) 



In the region of interest, where R — t- Rh and therefore 
/ — !■ 0, the kinetic term in Sin dominates over that in 
Sout while the gradient term in Sin is subdominant to 
that in 5*0114. This yields the approximate action 



S K 2ti dt 



Rh 



drr 



drr'^f{dr^f 



Rh 



(14) 

Substituting the expansion for the scalar field <& produces 



dt 



^" , 2Sfe,fe'afc(i)a/c'(t)ufe(r)?ifc'(r) 



~2n I drr 
'0 



/•OO 

27r / drr^f Vafe(i)afc'(t)ufc(r)Mfc,(r) 



Now let 



and 



Mkk' ^4,^1 drr Uk{r)uk'{r) 

'0 



Nkk'^-^TT drr^f{r)u'k{r)u't,,{r), 
JRh 



which allows us to rewrite ( 15 ) as 

s^Y^ht 

1, "J 



k.k' 



—dkMkk'dk' — -akNkk'ik' 



The action ( 18 1 corresponds to the Hamiltonian 



(15) 



(16) 



(17) 



(18) 



(19) 
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where 



n 



5L 

da 



(20) 



Since M and N are Hermitian matrices, it is possible to 
do a principle axis transformation to diagonalize them 
simultaneously. The exact form of them, however, is not 
important here. The Schrodinger Equation for a single 
mode will then be 



- — / — 



1 



Kb^ 



dt - 



(21) 



where m and K are the eigenvalues of M and N, respec- 
tively, and h is the eigenmode. 

This, however, only represents the Hamiltonian of the 
radiation, not the total system. The Hamiltonian of the 
system can be written as 



ill 



Ref.f 



Htotal — Hu 

showed that 

Hlaii = imu? + /(2vr/ii?")', 



(22) 



(23) 



where ^, — — 2TxaGRH and a is the tension in 

the wall. As i? — Rh, / — > 0, and the first term on the 
right hand side will dominate (since 11 ^ /~^''^), yielding 



wall 



(24) 



where we have chosen the negative solution because the 
shell is collapsing. The Hamiltonian of the radiation is 



/ 



K 



2to " 2 

The Hamiltonian of the entire system is therefore 



total 



2m' 



H 



(25) 



(26) 



Note that the momentum operators have different in- 
dices, where H^j = —id/dR is the momentum operator 
conjugate to the shell's position and Hf, = —id/db is the 
momentum operator conjugate to the eigeinmode b. The 
wavefunction of the entire system is therefore a function 
of b, R, and t, which we will write as 



* = *(6,i?,i). 



(27) 



Using (26), the Functional Schrodinger equation can be 



written as 



9* 

dR ' 2m 



9^- 



(28) 



Solving this requires knowing the classical equation of 
motion for the distance of the domain wall, which we 
will take as approximately 



R' 



(29) 



where the negative sign was chosen because the shell is 
collapsing. Since R(t) is only a function of time, we can 
rewrite ([28| as 



..19* / K 



R dt 2m 962 2 



94- 



According to (291, however, f /R = — 1, therefore 

/ 92* X,,^ „ 9* 
- 1 fe2^ ^ 2i — . 

2 at 



2m 952 



(30) 



(31) 



Rewriting this in the standard simple harmonic oscillator 
form yields 



2m 9fe2 



'"■ 2.2 

— uj b 

2 



-qi^b, r/) = i 



9*(&,77) 
di] 



where 



and 



.2 _ K _cji 



mf f 



(32) 



(33) 



(34) 



Here we have chosen rj{t 0) = 0. 

At early times, the initial vacuum state for the wave- 
function %p is that of a simple harmonic oscillator given 
by 



ip{b,r) = 0) 



mojQ 



1/4 



(35) 



At later times, the exact solution for the wavefunction is 
(see Ref.[T3j) 



ij{b,r]) 



1(77) 



1/4 



exp 



(36) 

where prj denotes the derivative of p{ri) with respect to 
77, and p{ri) is the solution to the equation 



d^P , 2/ ^ 1 



with initial conditions 



p{0) = 



P.(0) = ||o=0. 



The phase a is given by 



a (77) = - 



1 /•" dif 

2 Jo pHv')' 



(37) 

(38) 
(39) 

(40) 
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IV. OCCUPATION NUMBER OF THE 
RADIATION 

Consider an observer with detectors that are designed 
to register particles for the scalar field (p at early times. 
At late times, the observer will interpret each mode b 
in terms of the simple harmonic oscillator states, with 
final frequency w. The number of quanta in eigenmode 
b can be found by decomposing the wavefunction ( 36 ) 



into the simple harmonic oscillator states and evaluating 
the occupation number. The wavefunction -0 written in 
terms of the simple harmonic basis, {(pn}, at i = is 
given by 



where 



C„(t) = / dbipl{b)ij{b,t), 



(41) 



(42) 



which is the overlap of a Gaussian with the simple har- 
monic basis functions. The occupation number at eigen- 
frequency uj is given by the expectation value 



(43) 



After substitution, we find that the occupation number in 
the eigenmode b is given by (see Appendix B in Ref.|14j') 



(44) 



We will find that the plot is very similar to a typical 
Planck distribution. Therefore, we can find the tempera- 
ture of the radiation by treating the occupation number 
of each eigenmode b as if it followed the Plank distribu- 
tion. 



1 



(45) 



For the total system, the plot of occupation number 
as a function of frequency is shown in Fig.Q. Notice 
how similar this is to a Planck distribution. In fact, as 
time goes on, the distribution becomes more and more 
Planck-like as the oscillations in the curve dampen out 
and the tj = occupation number gets larger. If one 
plots In (1 + 1/A^) as a function of (I;, the slope will yield 

The plot of ln{\ + 1/iV) as a function of frequency is 
shown in Fig. ([2]). The straight lines represent best-fit 
lines whose slopes represent /3 at the chosen time. It 
should be noted at this point that /3 is in terms of scaled 
time T] instead of the asymptotic observer time t. There- 
fore, we must rescale /3 back to the observer's time, which 
is achieved by 



/3 



(46) 



While in scaled time f] the /3 actually diverges, in the 
asymptotic observer's time t the /3 actually approaches a 
constant, as shown in Fig.([3|. 

Finally, the temperature as a function of the observer's 
time can be found by inverting /3, as shown in Fig.Q. 



T(t) = 



(47) 
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FIG. 1: Occupation number of the radiation as a function of 
frequency. From bottom to top, the three lines are for t=13, 
t=14, and t=15, respectively. 




500 1000 1500 2000 2500 3000 
FIG. 2: The best-fit slopes of these curves represent /3. 

Notice how, from Fig.Q, there is a short time where 
the temperature increases before it starts to decrease 
again. This is an effect which is much more prominent 
when the domain wall starts just inside the cosmological 
horizon, and it is suppressed the smaller the cosmological 
constant is. Notice also that the temperature approaches 
a constant at late times. According to [T7], the Hawking 
temperature of a de Sitter Schwarzschild black hole is 



H 



1 

8nM 



AM 
~2^' 



(48) 
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FIG. 3: /3 as a function of time for a collapsing de Sitter 
Schwarzschild domain wall. 
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FIG. 5: Temperatme as a function of mass for a de Sitter 
Schwarzschild domain wall. 
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FIG. 4: Temperature as a function of time for a collapsing de 
Sitter Schwarzschild domain wall. 



From this we find that the ratio of the late time temper- 
ature in Fig.Q to the accepted Hawking temperature 
is Tschwh/TH = 1.02. We will now determine how the 
late time temperature scales with the mass of the domain 
wall. The plot of the domain wall's late time tempera- 
ture as a function of mass, along with the best-fit curve, 
is shown in Fig.([5]). The equation of the best-fit curve is 



SchwA 



0.0397304 

M 



0.0138729M. 



(49) 



Notice that the scaling with mass shows very good agree- 
ment with the expected result. 



V. ENTROPY 

The thermodynamic definition of entropy in terms of 
temperature is 



5 = 



dQ 
T ■ 



(50) 



Since changing the energy of the domain wall is the same 
as changing the mass, this can also be written as 



dM 



(51) 



Therefore, if we can determine how the temperature of 
the domain wall depends on the mass, we will be able to 
find an expression for the entropy in terms of the tem- 
perature. And since we know the time evolution of the 
temperature, we can determine the time evolution of the 
entropy as well. 



We insert Eq.(49) into Eq.(51|, and after integrating 



we find that the entropy as a function of temperature is 



S 



7 



12iy3(T- z/) 



372 i/A 



6YA 



+ 2^(-7-^+7^A)/n(7+^), 



(52) 



where 7 = 0.0397304 and 1^ = -0.0138729. After normal- 
izing the entropy such that the initial entropy is zero, we 
find that the entropy approaches a constant, as shown 
in Fig.Q. We also wish to point out here that, though 
the effect is small, the entropy does in fact decrease for a 
short time before increasing. The time dependent analy- 
sis of the 3-1-1 BTZ domain wall also exhibited a decrease 
with entropy over time, albeit in that case the effect was 
ever-present instead of temporary as it is here. This im- 
plies that the probable cause of the decreasing entropy is 
the cosmological constant and not the difference in topol- 
ogy. The decreasing entropy also means that for this 
particular choice of parameters, the domain wall will not 
collapse spontaneously, but rather would require an en- 
ergy input to get it started. The explanation of this effect 
is quite intuitive, as the collapse is at odds with the ex- 
pansion of space, and if the cosmological constant is large 
enough, then collapse will be prevented. We discovered 
that this feature is a function of the particular parameters 
that we chose. The amount of decrease is amplified when 
the domain wall starts close to the cosmological horizon. 
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FIG. 6: Entropy as a function of time for a de Sitter 
Schwarzschild domain wall. 

and there is no longer a decrease in entropy when the cos- 
mological constant is sufficiently small. Given how small 
the cosmological constant is by current measurements, it 
is unlikely that the collapse of a de Sitter Schwarzschild 
domain wall would be forbidden in our universe. 



VI. CONCLUSION 

We investigated the time evolution of the temperature 
and entropy of a de Sitter Schwarzschild domain wall 
(representing a spherically symmetric domain wall with 



a cosmological constant). This was done by coupling a 
scalar field to the background of the domain wall and 
evaluating the occupation number. From the occupation 
number we were able to determine /3 and therefore the 
temperature. We found that the temperature of the de 
Sitter Schwarzschild domain wall increased for a short 
time before then decreasing and approaching a constant. 
We found that this effect was more prominent when the 
initial position of the domain wall was close to the cos- 
mological horizon. The effect is also suppressed when the 
cosmological constant is small. The late time tempera- 
ture exhibited very good agreement with the Hawking 
temperature of a static black hole. It also produced scal- 
ing with mass that matched the literature. The entropy, 
however, was interesting because it exhibited qualities 
that were seen both in the Schwarzschild domain wall 
and in the 3-1-1 BTZ domain wall. Namely, overall the 
entropy increased and approached a constant, but there 
was also a short time in the beginning of the collapse 
when the entropy decreased by a small amount. Com- 
parison with the results for the Schwarzschild and BTZ 
domain walls leads us to posit that the cause of the de- 
crease in entropy was the cosmological constant and not 
the difference in topology. 
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